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We provide analytic solutions for the power spectrum and bispectrum of curvature fluctuations 
produced by a step feature in the inflaton potential, valid in the limit that the step is short and 
sharp. In this limit, the bispectrum is strongly scale dependent and its effective non-linearity 
attains a large oscillatory amplitude. The perturbations to the curvature power spectrum, on the 
other hand, remain a small component on top of the usual spectrum of fluctuations generated by 
slow roll. We utilize our analytic solutions to assess the observability of the predicted non-Gaussian 
signatures and show that, if present, only very sharp steps on scales larger than ~ 2 Gpc are likely 
to be able to be detected by Planck. Such features are not only consistent with WMAP7 data, 
but can also improve its likelihood by 2AlnL ~ 12 for two extra parameters, the step location and 
height. If this improvement were due to a slow roll violating step as considered here, a bispectrum 
or corresponding polarization power spectrum detection would provide definitive checks as to its 
primordial origin. 



I. INTRODUCTION 

In this paper we employ the generalized slow roll 
(GSR) approach [IHJ to consider the bispectrum of pri- 
mordial fluctuations generated by a transient violation 
of slow roll due to a step feature in the inflaton poten- 
tial. While earlier works in the field have considered the 
non-Gaussianity from these models, here we explore a 
hitherto unexplored region of the parameter space where 
the step in the potential becomes sharp. This region of 
parameter space is near the so-called "decoupling" limit 
where the perturbations in the power spectrum due to 
the step remain small while higher order TV-point func- 
tions are parametrically enhanced [SJ |5] . 

Features in the inflationary potential have a long his- 
tory [7]. They came into vogue as a possible explana- 
tion for the apparent low multipole glitch at ^ ^ 20 — 40 
in the angular spectrum of the cosmic microwave back- 
ground (CMB) radiation |5HTT]. If such a feature arises 
as a result of a glitch in the primordial potential then 
corresponding features should arise in the electric-type 
polarization pattern of the CMB [11]. More generally, 
failure to detect consistent patterns in the polarization 
may even be considered evidence against a single adia- 
batic degree of freedom as the source of the primordial 
perturbations [4]. 

Beyond two-point statistics such as the angular tem- 
perature and polarization power spectra, it is well known 
that such a glitch in the potential leads to higher A^-point 
functions which are enhanced relative to those generated 
by the otherwise smooth background [12nl6j . However, 
as was shown in our earlier work |17) . the bispectrum, 
or three-point function produced by a step feature that 
best fits the glitch in the WMAP data (H HI] falls short 
of detectability by a very wide margin. 

In this work we investigate the space of parameters 



that correspond to the step becoming short and sharp. 
As we will see, in this limit, the perturbations in the 
power spectrum remain small while the reduced bispec- 
trum or effective non-linearity, /nlI^) becomes large. 
Despite this seeming largeness, the bispectrum remains 
difficult to detect due to its oscillatory form. Nonetheless, 
for a sufficiently sharp step at sufficiently large scales, an 
observably large non-Gaussianity is compatible with cur- 
rent power spectrum constraints. 

This paper is organized as follows: In Sjlljwe derive ap- 
proximate solutions to the evolution of the inflaton on an 
inflationary potential that undergoes a sharp downward 
step. These solutions yield analytic expressions for the 



perturbation to the curvature power spectrum in \ III 



In iV\f we review the formalism used to calculate the 



bispectrum from which we obtain approximate analytic 
solutions. The constraints that the WMAP 7 year power 
spectrum data place on the sharp step region of parame- 
ter space in ^jVj inform the prospects for the detection of 
the bispectrum in a cosmic variance limited CMB exper- 
iment in §VI| We summarize our findings in j jVII[ 

Some related but separate studies are presented in the 
Appendixes. In Appendix |A] we investigate the real 
space analogs of the power spectrum and bispectrum of 
curvature fluctuations generated by a sharp-step feature. 
In Appendix [B] we calculate the leading order slow roll 
corrections to our analytic solutions. Finally in Appendix 
[C| we demonstrate that our analytic solutions are approx- 
imately separable. 

Throughout the paper we work in units where the re- 
duced Planck mass Afpi — (SttG)^^/'^ = 1 as well as 



II. BACKGROUND SOLUTION 

We consider a homogeneous scalar field minimally cou- 
pled to gravity with a potential of the form 



1/(0) = ^0(0) 



l + cF 



(1) 



where the potential Vq supports slow roll inflation in the 
limit where c — 0, and F{x) is a function which tran- 
sitions from —1 to +1 as its argument passes through 
a; = from above with a characteristic width Ax — 1. 
Thus the potential characterizes a step of height 2c and 
half-width d. When an explicit choice is required, we will 
take 



F{x) = — tanh(2;), 



(2) 



as an example in the following sections. 

The scalar field obeys the usual equation of motion on 
the background 



dV 



= 0, 



(3) 



where H = d/a is the Hubble parameter, and an overdot 
represents a derivative with respect to cosmic time. 

In the limit in which the step is short, c <C 1, its 
presence in the potential will have only a small effect 
on the evolution of the inflaton. Specifically, as long as 
the change in the potential as the field rolls across the 
step is small compared with the kinetic energy of the 
background field, then we can think of the effect of the 
step as a perturbation on the evolution of on the po- 
tential Vq. As the field crosses the step a potential energy 
Ay cV is converted into kinetic energy. Comparing 
this with the kinetic energy of the inflaton rolling on the 
background we find 



2AV 



^2 



6- 



where we have introduced the slow roll parameter 

12 



en 



Whenever 



H 

c < 



(4) 



(5) 



(6) 



we can treat the effect of the step as a small perturbation 
on the background described by Vq. 

With this limit in mind, we look to solve Eq. ([s]) itcr- 
atively as follows: taking 6c <C en, we split the field into 
a piece zeroth order in c and pieces which we take to be 
higher order in c 



+ 



en = €a + ei 



(7) 



In this expression, 0o characterizes the behavior of the 
field in the absence of the step, i.e. on the potential Vq, 



while 01 is the perturbation which is taken to be linear 
in c, and '. . .' denotes terms higher order in an expansion 
in 



C = 



6c 



(8) 



While the condition in Eq. ([6|) requires that 0i remains 
small, it places no restriction on higher derivatives. In 
particular, we will see that by making the step sharp 
we can make the acceleration and jerk of the field 0i 
arbitrarily large. 

It proves useful to write the equation of motion for the 
perturbation, 0i , using the value of the background field 
00 as the independent variable. The condition 6c ^ eg 
requires that (f)o =/= and 4>oit) evolve monotonically in 
the vicinity of the step which ensures that 0o is a suitable 
time variable. Given that on the background potential 
Vq, the slow roll approximation yields eo ~ {V /V)'^/2 
and assuming that eq ~ const, (and so not ^ 
we obtain the equation of motion for 0i as 



d 
d0o 



-F 



dF 



(9) 



where the subscript q denotes quantities to be evaluated 
on the unperturbed (c — 0) background and we have 
dropped some terms that are suppressed by slow roll pa- 
rameters evaluated on the unperturbed background. 

For small sharp steps, the dF/d(j)Q term in Eq. Q dom- 
inates and 



dcpn 



const. 



(10) 



The exponential in the integrand is slowly varying across 
the width of the step and may be evaluated at 0/ and 
taken out of the integral. Furthermore the integration 
constant is set by the boundary condition that d(j)i /(i0o — 
before crossing the step and so 



l-F 



(11) 



how- 



In principle, we could further integrate to find 0i; 
ever, this is not necessary for this work. 

The final step is to reference this solution to conformal 
time rather than the background field. Since the pertur- 
bation 01 is only important for a short period of time, 
the variation of the unperturbed (c = 0) slow roll param- 
eters is negligible. In this approximation, the evolution 
of the background in conformal time is very simple 



d0o 



= ±v2eo ~ const.. 



din 77 

t>o-<l>f = ±\/2eo ln(?//?y/). 



(12) 



where the sign should be chosen depending on the direc- 
tion the field is rolling - i.e. whether one is in a large, or 
small field inflationary model and one is rolling toward 



2 



or away from the origin. In this work, we have large field 
inflationary models in mind, e.g. Vq = ^m'^cp'^ and will 
take the positive sign. 

With these relations, we can explicitly solve for the first 
order changes to the slow roll parameters as a function 
of conformal time. The slow roll parameters are defined 
by Eq. ^ and 



(j)H 



So = 



(13) 



Using the results above, we divide the slow roll parame- 
ters into a background part and a part that is linear in 
the perturbation, rjH = Vo + Vi ^-nd S2 = ^2,0 + '^2,11 



C 



1 - F 



In 



(14) 



Here and below primes denote derivatives with respect 
to In 77. Notice that ei remains a tiny correction to the 
already small ep, of order 0{C), as expected. On the 
other hand, in principle rji and ^2,1 can be arbitrarily 
large due to the presence of the 1/d factors from the 
derivatives. In what follows we will see that these give 
rise to a large and strongly scale dependent bispectrum. 

In deriving these results, we have assumed that the 
perturbation series converges. At first, it may seem that 
this requirement would impose a restriction on the width 
of the step, d, due to successively expanding the dV/d(f) 
term in Eq. ^ at each order in perturbation theory. 
However, note that by iteration, these terms are always 
multiplied by the previous order in the perturbation se- 
ries 



\ d(t)Q I 



Cd"F 



-le 



(15) 



Around the step, the derivatives of F carry factors of d^" 
for a A(f)Q interval of order d. Thus in this neighborhood, 
the integrals over c/jq give contributions of order d, (pi 
is of order (Cd), and (/)„ is of order (Cd)". With these 
relations dcpn/dcpo is order C" there and is independent 
of d as d — )■ 0. The corrections to the slow roll parameter 
are then also strongly convergent as long as C <C 1. 

Another way of seeing why the series is convergent as 
we take d — >■ is to realize that these perturbative correc- 
tions amount to small shifts in when the inflaton crosses 
the sharp feature and hence unobservable changes in the 
temporal location of features in the slow roll parameters. 

In fact for the purpose of computing the shapes of these 
features, we can always define the position of the step in 
some arbitrary choice of the zero point in In 77 or efolds to 
be the same as that of the background to all orders. The 
background has an approximate time translation sym- 
metry (broken only weakly by ^ 0) and so there is no 
dynamical consequence to this choice. As d — )■ the step 
becomes sharp and changes to the potential are no longer 



well approximated by a Taylor series for finite perturba- 
tions in field value. But making use of the approximate 
time translation invariance of the background, we have 
defined these away. It then follows that the linear expan- 
sion in C remains valid in the limit d — > 0. We shall see 
in the next sections that nonlinear scalings of this type 
can be used to extend the analytic solutions to the power 
spectrum and bispectrum out to C ~ 1. 



III. CURVATURE POWER SPECTRUM 



In §111 A[ we use the background solution of the last 
section to derive an analytic form for the curvature power 
spectrum of the step potential. We test this analytic form 



against numerical calculations in § HI B and show that it 
can be extended to order unity features with a nonlinear 
rescaling of parameters. 



A. Analytic power spectrum solution 

In the generalized slow roll approximation, the power 
spectrum of curvature fluctuations = k^P-ji /2tt'^ is 
given to leading order in the source function by (2j 



InA^(fc) =G(lnr/,)+ / 

Jr), 



drj 



I^(fcr/)G'(lnr/), (16) 



where the source function 
G(ln77) 

with 



21n/+-(ln/)', 



/ = 



H 



{aHr]). 



(17) 



(18) 



Recall that primes denote derivatives with respect to In 77. 
To leading order in the slow roll approximation, note that 
is simply the power spectrum. The window function, 
W{x), is given by 



W(x) 



3sin(2a;) 3cos(2a::) 3 sin(2x) 



2x3 



2x 



(19) 



To linear order in the step height c or kinetic energy 
perturbation C, the power spectrum source is given by 
G' = G'o + G[ with 



G'i(ln77) = -(6ei- 3771+^2,1), 



(20) 



where Go(ln77) is the source function in the absence of 
the step. Compared to rji and 62^1, ci is negligible and 



F' 



(21) 



3 



We can now evaluate the power spectrum. Inserting 
Eq. ( 21 ) into Eq. ( 16 ), integrating by parts and dropping 



the neghgible boundary terms, we obtain 



C 



Vf 



F'. 
(22) 



In the d — Hmit, F' becomes 2 times a delta function 
at Inr^j^, given our convention that the step height is 2c 
and so in this limit 

hm In A^(fc) = In A^^„(/c) + ^W'{krjj). (23) 
Notice that since W'{x) — dW/dlnx, 
W'{x) = 



9 \ f 9 \ sin2a; , , 

3+^ C0S2.+ 15-^ — , (24) 



and in this limit there are oscillations of amplitude C in 
the power spectrum out to A: — >■ oo. This behavior is 
reminiscent of the Fourier transform of a function with a 
sharp feature in real space. In Appendix [Xj we shall see 
that this intuition is borne out by an analysis of the real 
space correlation function and that the high k oscillations 
correspond to a feature at r = 2r]f. 

As we move away from the d — >■ limit, the function F' 
has finite width of order Arj/rif ^ d/y/eo whereas the os- 
cillatory features in W' are of order Aij ~ 1/fc. Therefore 
the delta function approximation holds for krjf ^/eo/d. 
For larger values, we expect the integral to be smaller due 
to integration over the oscillations in W' . In other words, 
the 0? — >■ results should be multiplied by some damping 
factor that depends on krjfd/^JtQ. 

In order to explicitly evaluate this damping function, 
we now assume the form for F given in Eq. ([2]). The 
required integrals then take the form 



/ = 



dlnrj f{kri) sech^ 



V, 



(25) 



where f{krj) is a polynomial times a sine or cosine and we 
take r]^ to be a time well after the inflaton has crossed 
the feature. In the limit that d — >■ 0, the polynomial 
varies slowly while sech^ is non-zero, and consequently we 
can replace it by its value when sech^ is peaked, namely, 
[kr])^ — {krif)^. We are then left to evaluate integrals 
such as 



/ = 



din 77 sin(2fc77)sech^ 



■In 



(26) 



as well as the same thing with sine replaced by cosine. 
Working in the limit ry* and changing variable to 



In(^), the integral can be written 



dy sin 



2krif exp — 



d 



sech^y. 



(27) 




FIG. 1. Initial curvature power spectrum evaluated numerically 
(in black solid curves) compared with using the approximation in 
Eq. ^32\ (in dashed red curves) for a small amplitude step c = 
10~^ <SC eo/6 (or C <g; 1), where the approximation is expected to 
work, and a step width, d = 10"'^, or ft: 4.3. The lower panel 
shows the difference between the two curves, numerical - analytic. 
The remaining parameters used in both models are shown in Table 

m 



The integrand above only has support for y g (—1,1) 
and, in the limit d — > 0, the phase of the sine varies 
slowly, and thus we can expand. 



dy sin 



2k7Jf 1 



sech y. 



(28) 



The phase error of the expression in brackets is 
O^kjjf {yd/ ^/2eo)^) and, given the support, the approxi- 
mation is valid for krjf ^ 2neo/d'^. But, since d/i/2io <C 
1, this only breaks down for values of kr]f that have many 
oscillations in the region (—1, 1) and are hence far into 
the region where the integral is already negligible. The 
integral in Eq. ( pS] ) can then be performed by integrat- 
ing around the square contour in the complex plane, with 
vertices at (i?, 0), (i?, itt), {—R^m) and (— i?, 0) and tak- 
ing the limit i? — > c», obtaining. 



/«2^p(^)sin(2fc,,), 
and similarly for the terms with cos{2krjf). Here 



nd 



(29) 



(30) 



determines the value of krjf at which damping starts to 
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1 10 100 1 10 100 
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FIG. 2. Same as Fig. [Tlbut using a smaller value of the step width. 
Here we take d = 10^ , or z j 43. 



become important and 



sinhy 



(31) 



is the damping envelope for the tanh step profile. 
Generically then, the power spectrum takes the form 

In A^(fc) = In A^,o(fc) + (^^) W'ikrjf). (32) 

Here, the oscillating window function W' is modulated 
by the decaying envelope T> which is set by the details 
of the step. As a; — >■ 0, W'{x) — >■ 0, and no spurious su- 
perhorizon contributions during inflation are generated. 
Interestingly, the amplitude of the perturbation is inde- 
pendent of the width of the step, while the range (in k 
space) in which the perturbation persists is independent 
of the height of the step. This analytic solution was de- 
rived in a different way by Stewart [TH] who also noted 
that different functions F simply change the damping en- 
velope. Similar conclusions were also reached in connec- 
tion with steps arising from features due to duality cas- 
cades in brane inflation [19j . In Appendix [A] we address 
the paradox that even if the feature scale ?7/ is greater 
than the current horizon, the power spectrum can retain 
oscillations to arbitrarily high k. 



B. Numerical comparison and nonlinear scaling 

In Figs. [T] - [3] we compare this approximation to the 
exact solution. We choose a model with Vq = m^0^/2. 



FIG. 3. Initial curvature power spectrum evaluated numerically 
(solid black curves) compared with the approximation in Eq. | |32| l 
(dashed red curves) for a large amplitude c = eo/6 (or Co = 1, 
C = 2/3), where the linear expansion is violated by order unity 
terms. Using the leading order expressions for the step location, 
oscillation amplitude and damping r]f,C = Co and (upper panel) 
leads to phase and amplitude errors that are largely corrected by 
non-linearly rescaling these parameters using Eqs. | |33| l and | |35| l 
(lower panel). 



Parameter 


Value 


m 


7.126 X 10"® 


</>/ 


14.668 


eo 


0.00925 




2.231 




0.1114 


e 


1.0401 


r 


0.084 



TABLE I. Fiducial model parameters for a step at (/>/ on a 
potential Vb = m^0^/2 with physical baryon and cold dark 
matter densities Q,bh? , flch^, 100 times the angular size of the 
sound horizon 6, and the optical depth to reionization r. 



and show curves for d = 0.01 in Fig. [T]and d — 0.001 in 
Fig.[2| to illustrate the damping behavior. We summarize 
the other parameters used in Table III 

For c = 10"^ < eo/6 (see Figs, [l] and [2]), the agree- 
ment is excellent both in the amplitude and phase of 
the results. With c = 10~^ the inflaton crosses 0/ at 
7] = 1454.6 Mpc, which agrees to 0.1% with 77/ = 1456.1 
Mpc, defined as the epoch when the inflaton crosses the 
feature for no step c = 0. In the lower panels of Figs. 



5 



[T] and [2] we plot the difference between our analytic ap- 
proximation in Eq. ( 32 1 and a numerical evaluation of the 



spectrum. The error on our approximation here is at the 
level of 5%, and is contained in a component that is out 
of phase by 7r/2 with our approximation. As we point out 
in Appendix [B| this error can be attributed to slow roll 
corrections to the modefunctions which are not captured 
by our leading order approximation at Eq. ( 16 ). Further- 



more, the spectra also disagree on the percent level on 
scales far from the k ~ ^/Vf- This disagreement is due 
to our analytic approximation of Eq. (16 1 by Eq. (32) 



and can be corrected in a straightforward manner by in- 
cluding subleading terms in the power spectrum source, 
G', at Eq. 

For c = e^j6 (see Fig. |3j upper panel), we see both an 
error in the amplitude, especially at the first few oscilla- 
tions, as well as a difference in phase. As discussed in the 
previous section, the phase error simply represents a dif- 
ference between when the inflaton crosses the feature at 
(f>f relative to the c — model. Here the inflaton crosses 
the feature sXrj — 1271.7 Mpc. If we instead define 77/ in 
the analytic expression Eq. (l32|) to be 



(33) 



T]f = ri{(l)f), 



we eliminate the phase error. We can also substantially 
improve the fit by rescaling of the value of eo to include 
the first order perturbation as the field crosses the step, 



<^HiVf) ~ eo + 3c. 



(34) 



Note that the first order expression suffices since its ac- 
curacy is guaranteed by energy conservation so long as 
c ^ 1. In particular we can define the non-linear ampli- 
tude and damping as 



C = 



6c 
eo + 3c ' 
V2eo + 6c 
7rd 



(35) 



We demonstrate these non-linear corrections with the 
general form of Eq. ( 32 ) in the lower panel of Fig. [3] The 
parameter C in the lower panel is related to the original 
linearized version Cq = 6c/eo in the upper panel by 



C = 



l+Co/2' 



(36) 



We employ these definitions of 77/, C and Xd from this 
point forward. In Table |llj we compile a guide to the 
step parameter notation used in the rest of this paper. 



IV. CURVATURE BISPECTRUM 



In \ IV A| we briefly describe the method we use to cal- 
culate the bispectrum. For more detail, we refer the 
reader to [T71 [201 HI]- We use this method to derive 



an analytic approximation in S IV B In S IV C we discuss 



Parameter Definition 



Eq. 



F Step shape function 

(j>f Potential step position 

rjf Step crossing time 

c Potential step height 

C Kinetic energy perturbation 

Co Linearized C 

Ac Angular power spectrum amplitude 

Ab Angular bispectrum amplitude 

d Potential step width 

Xd Dimensionless damping scale 

Id Angular damping multipole 

2? Damping function 



TABLE II. Step parameters and their defining equations. 



limiting cases and consistency checks for the approxima- 



tion and in ^VD we compare it with direct numerical 



computation of the exact solution. 



A. The action and the in-in formalism 

In this paper to move beyond the linear equations of 
motion, we work in comoving gauge in which the time 
slicing is chosen so that the inflaton is unperturbed. In 
this gauge the scalar degree of freedom is the comov- 
ing curvature perturbation, while the remaining physical 
metric degrees of freedom are the two polarizations of 
the transverse traceless tensor perturbation. Tensor per- 
turbations are not enhanced due to the presence of the 
step ^22, and thus we neglect them here. We make use 
of the interaction picture, where we choose our basis so 
that the flelds diagonalize the Hamiltonian arising from 
the quadratic action. 



^2 



(37) 



In this picture the Hamiltonian arising from cubic and 
higher order terms in the Lagrangian, which we shall refer 
to as the interaction Hamiltonian, Hj, evolve the states. 
At leading order, the tree-level bispectrum at time is 
then given by [20] 



(7^kl(^*)7^k.(^*)7^k3(^*)> 



23? 



(7^kl {t, )7^k, (t* )7^k3 {u)Hi (t)) 



(38) 



At leading order in fluctuations, the interaction Hamil- 
tonian in Eq. ( 38 ) arises from the cubic action. In this 



work we are interested in potentials in which the inflaton 
undergoes a sharp transient acceleration but inflation is 
not interrupted. Consequently, e// <C 1 everywhere and. 



6 



to a very good approximation, the cubic action is given 
by 



dtcfx 



ahnien -m)TZ'''JZ 



dt 



(39) 



Note that, after switching to conformal time, 77 = 
J^""'' dt' / a{t) (defined to be positive during inflation), 
and performing some integration by parts, the action in 
Eq. ( 39 ) can be written as 



and thus (for this interaction), 

/:3 = 2(77H-eH)7^/:2. 



i&Rf] (40) 



(41) 



Naively, in order that we trust our perturbative evalua- 
tion of correlation functions using the interaction picture, 
we might demand that £3 < C2 so that the evolution of 
the operators is well described by the equations of mo- 
tion arising from the quadratic action. Thus we require 
at least 



2(r]H - eH)'R < 1. 



(42) 



We shall see that in the context of the step model, this 
constraint places a lower limit on the step width d that 
we can consider using perturbative techniques. 

In conformal time, the relevant interaction Hamilto- 



nian that follows from the action in Eq. ( 39 ) is given by 
d^qa d^qb d^qc 



(27r)3 (27r)3 (27r)3 



{27ryS^{(la + qb + qc) 



3ry2 



(e^f - 7;h)' 7^q„7^q,7^qJ (43) 



dr\ \ Srj 



(ei^-7/J^)(7^q„7^q,7^qJ' 



where recall ' = d/dlnrj. In this expression, the fields 
TZ are interaction picture fields chosen to diagonalize 
the Hamiltonian derived from the quadratic action in 



Eq. (37) 



We define the bispectrum through 
where 

r^a'eHieH-VHnnlTZlTZlJ 



(44) 



7-1=7], 



(45) 



and 3ff denotes the real part. We choose 77* to be a time 
well after the inflaton has crossed the feature, and ei = 
i]i =(^2,1 =0. Consequently, the second term in Eq. (|45|) 



is tiny, and we neglect it in what follows. The derivation 
of Eq. (451 is described in detail in [T7] , and we refer the 
reader to this work for details. 



B. Analytic bispectrum solutions 

As described in [17], the GSR approach proceeds by 
iteratively correcting the evolution of the mode func- 
tion for the effect of deviations from de Sitter space. 
As we will see in this section, in the limit that d 
the modefunctions, TZk, are well treated as unperturbed, 
and the non-Gaussian features are well described by in- 
tegrals involving only the unperturbed modefunctions. 
We have already seen that, even in the limit d — ?> 0, the 
power spectrum, which may be thought of as the square 
modulus of the modefunctions, gains a correction that is 
O^c/en)- Thus at linear order, the corrections involving 
the perturbed mode functions will be small compared to 
those arising from the perturbed slow roll parameters. 

There is an additional weak constraint on d from 
Eq. (42) due to the perturbative expansion. The per- 
turbations in the curvature themselves do not become 
large as the field crosses the feature, and thus we may 
approximate TZ ^ 10^^. Since en remains small, we can 
neglect it and focus on rjn- Near the step, from above, 
we know that 



c 1 
e^d' 



(46) 



Thus, with eg ~ 0.01, this implies the limit on the ratio 



< 10* 



(47) 



In the limit of infinitesimal width d — >■ with finite 
height, while the field fluctuations remain small, they are 
no longer well characterized by their tree level iV-point 
functions. On the other hand, this is not a particularly 
restrictive constraint and when considering the impact 
on the primary CMB anisotropy is equivalent to taking 
d — > 0. Nevertheless, it does cure the weak logarithmic 
divergence in the real space correlation function discussed 
in Appendix [K\ 

In the GSR approach, as described in detail in [IT], we 
can approximate Eq. mm for the bispectrum as 



Bn{ki,k2,ki 



(27r)4 ^n{ki)^n{k2)l^n{kz) 



2^3 7^3 z^3 
""l "-2 ""3 



- Io(K)kik2k^ ~ h{K)Y,kikj 



+l2{K)K{kl + kl + kl)], (48) 

where = /ci 4- A;2 -|- ^3 is the perimeter of the triangle in 
momentum space. The principle advantage of Eq. ( 48 ) is 
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that it involves integrals universal in K , 

(49) 

where 



Wo{x) 



xsmx. 



cos a;, 



W2{X) 



smx 



(50) 



and the bispectrum source is given by 

Gb = (^) . (51) 

The function A7j(fc) is given by the root of the power 
spectrum, ATi{k) — ^A^(fc). 

We have already obtained analytic solutions for the 
power spectrum A-iz{k) above, and it remains to evalu- 
ate the integrals Iq, Ii and l2- Following the analysis for 
the power spectrum, we split the bispectrum source into 
a background piece and feature piece, but here the back- 
ground piece generates a negligible bispectrum. Likewise 
the boundary term in Eq. ( 49 ) is negligible and the 



term in the source is suppressed by a factor of eo in its 
contribution. Integrating twice by parts, we then obtain 
from Eq. (Ill I 



= ^ r^(^^-l)f^) (52) 



We can bring this into the form of Eq. ( 22 1 with one more 



integration by parts that defines a new window function 

X^{K^) = ~ fnYw^iKfj). (53) 

Jo V \VfJ 



Using the approximation for the tanh step in Eq. ( 29 ) we 
obtain 

=1;^ (Sf ) 



where recall Xd was defined in Eq. ( 30 1 such that x^ ^ oo 
as d — > 0. Explicitly 



Xoix)=- 



{x'^ - 9x^ + 54) 
x"^ 

„4 0-7^2 



COS X 



2x^-27x^ + 54) sinx 



3 (x^ — 6) cosx [x^ 



Xi{x) = 2 ■ 3 

X'^ X'^ 

(x^ — 9) cosx (4a;^ — 9) sinx 



6) (x^ — 3) sinx 



X2{x) 



(55) 



All of these window functions vanish as X as X y 
0, which means that we do not generate any spurious 
superhorizon effects. Notice also that Xq diverges as x^ 
for large x while Xi diverges linearly and X3 approaches 



a constant in this limit. This implies that in the limit d 



the quantity in brackets in Eq. ( 48 ) increases without 
bound as fc — > 00. For small but finite d, we can estimate 
the location and height of the peak non-Gaussianity. In 
this limit the bispectrum is dominated by the quadratic 
term Iq, 



(56) 



whose envelope behaves as x'^ exp(— x) near its peak x 
3 or 



-f^'pcak = 6 



Xd 

Vf' 



(57) 



For a fixed step position, i.e. a fixed value of 77/, as c? — 0, 
the peak of the reduced bispectrum moves to larger and 
larger values of K. Simultaneously reducing the width 
of the feature and moving it to larger scales, i.e. increas- 
ing rjf, fixes the position of the peak of the reduced bis- 
pectrum, while increasing its amplitude proportional to 
the square of the ratio of step positions and increasing 
the frequency of its oscillations. Thus the sharper the 
feature, the larger the non-Gaussianity as anticipated in 

[B 113 [Ml- 

Indeed naively, it might seem that the non-Gaussianity 
strongly diverges with decreasing d due to the high k- 
modes. Here one must be careful in defining what is 
meant by a large non-Gaussianity. Note that the bis- 
pectrum itself, i.e. the amplitude of individual triangles 
in /c-space, does not diverge. It is only when written in 
terms of the reduced bispectrum or effective /nli that 
the apparent quadratic divergence appears 

lu h BTi{ki,k2,k3) 

^/NL(fcl,^2,fc3) =^-777^^-777^— , (58) 

5 Pn[ki)Pn[k2) +veTin. 

where "perm." refers to cyclic permutation of the indices. 
To see where this arises, we can rewrite the reduced bis- 
pectrum in terms of the quantity defined by |12[ I13j 



giki,k2,k3) _ klklk 



kik2k3 



4-BTiiki,k2,k3) 



(2^)4A 

,l^fc[+fc| + fc|\ Q.(^f^^^k2,k,), 
fcifc2fc3 / 5 



(59) 



where Ag is a constant taken to be Ag — for nearly 
scale invariant spectra. The extra factors of k account for 
the volume factors in k space or equivalently the number 
of fc-space triangles available. Divergence in the reduced 
bispectrum is usually associated with an equivalent di- 
vergence in the non-Gaussianity by assuming that the 
individual triangles add coherently. For oscillatory bis- 
pectra, counting triangles through k factors leads to a 
misleading sense of both the UV divergence and the ob- 
servability of the non-Gaussianity. 

A direct means of seeing this fallacy is to evaluate the 
non-Gaussianity in the real space three-point correlation 
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function. We show in Appendix that, similar to the 
power spectrum oscillations, the oscillating and diverging 
reduced bispectrum is associated with sharp features in 
the three-point correlation on scales r ^ rjf whose ampli- 
tude diverges no more than logarithmically. The oscilla- 
tions in the bispectrum prevent triangles from coherently 
adding up to a large non-Gaussianity in real space. Fur- 
thermore we have seen that there is a mild limitation 
on our derivation oi c/d < 10"^ from the validity of the 
perturbative expansion of the action which prevents even 
the logarithmic divergence from being manifest. 



C. Squeezed limit 

It is well known that the bispectrum of curvature fluc- 
tuations produced by a single scalar field during infla- 
tion satisfies a consistency relation which relates the bis- 
pectrum in the squeezed limit to the slope of the power 
spectrum [201 US] ■ Since we have analytic solutions for 
both the bispectrum and power spectrum, this provides 
us with a non-trivial test of our results. 

In the limit where one of the momenta is much smaller 
than the other two, ks ^ kL, the consistency relation 
implies for the reduced bispectrum 



6 



7/NL(fcs,fcL,fcL) 

5 



1 dlnA^ 



2 dink 



,g(fcs,fcL,fcL) 



(60) 



where the approximation holds for nearly scale invariant 
power spectra. Equation (48 1 then implies that, in this 
limit, 



din 



din A; 



(61) 



Now, from Eq. ( 32 1 , and dropping terms of order 
O(d/^2eo) - which amounts to ignoring the variation 
of the envelope function - we find 



din 



din fc 



n,{ki^)-l + ^V 



fcL 



kr]f\ 



W"{kT]f). (62) 



Comparing Eq. (62 I to Eq. (61) using Eq. (55), and ig- 



noring the slow roll suppressed terms, we find that in this 
approximation, the consistency relation is satisfied. 
A second non-trivial check of our results is that correc- 



tions to the squeezed limit of the bispectrum in Eq. ( 48 ) 
are proportional to (^s/^l)^; where fcg is the side that 
is being taken to zero. This is in accordance with the 
behavior of the so called "not so squeezed limit" [26] . 



D. Numerical comparison 

We begin with a comparison for a small step ampli- 
tude. In Fig. [4] we plot the bispectrum resulting from a 



2000 - 




FIG. 4. Equilateral bispectrum computed using the analytic ap- 
proximation in Eq. | |48| (red dashed) compared with a full numer- 
ical evaluation of the integral in Eq. l |45| for a small amplitude 
c = 10~^ (or C <K 1). The lower panel shows the difference be- 
tween the curves, expressed as a fraction of the envelope of the 
analytic approximation in Eq. K3t. 



step with height c = 10 ^ and width d = 0.001. We show 



the result of evaluating Eq. ( 45 ) numerically as well as 



the result of using the analytic approximation in Eq. ( 48 ) 



Overall, in both amplitude and frequency the approxima- 
tion is excellent. We show a more detailed comparison 
in the lower panel, where we divide the difference by the 
envelope of the analytic expectation 



9 ^koik) 



8 A 



3/2 



(63) 



where the amplitude of the powerspectrum is taken 
throughout this section to be As = 2 Ah x 10~^. The 
good overall agreement hides an oscillatory 10% residual 
error. 

Given that we are perturbing in C 0.6% here, this 
cannot be an O(C^) correction. In Appendix |b] we 
demonstrate that these corrections are largely due to slow 
roll corrections to the modefunctions. On superhorizon 
scales, these corrections lead to a phase shift in the grow- 
ing mode of the curvature perturbation, resulting in a 
correction on the 10% level to the bispectrum. In fact, 
as pointed out by |27j , leading order slow roll corrections 
generically lead to 10% rather than ©(eq) ^ 1% correc- 
tions as one might naively guess. Additionally, note that 
the error begins to increase for modes that are far into 
the damping window. These momenta have many oscil- 
lations during the period where the bispectrum source in 



Eq. ( 49 ) is non-zero and are thus in the region where our 
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FIG. 5. Equilateral bispectrum computed using the analytic ap- 
proximation in Eq. ( |48| (red dashed) compared with a full numer- 
ical evaluation of the integral in Eq. l |45[ l for a large amplitude 
c ft! eo/6 (or Co = 1,C = 2/3). Here the location, amplitude and 
damping scale of the feature have been nonlinearly scaled according 
to Eq. ^ and Eq. ( [35| . 
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FIG. 6. Fractional difference of the temperature power spectrum 
of a model with rif = 1.44 Gpc, c = 10~^ and d = 10~^, relative 
to a model with c = per unit power spectrum feature amplitude 
Ac (see Eq. |69| l . The case without gravitational lensing is shown 
in red dashed lines, and when gravitational lensing is taken into 
account is shown in black solid lines. In blac k d ot-dashed lines we 
show the expected envelope from Eq. |[70l. 



approximation in Eq. ( 29 ) is beginning to break down. 



Given that this is occurring in the region where the bis- 
pectrum is strongly damped, it will have little effect on 
our results. 

As we increase the step height, the leading order ap- 
proximation for r]f, C and Xd in Eq. (53) breaks down. 



Just as in the power spectrum, these errors in phase and 
amplitude are largely corrected by the nonlinear rescal- 
ings of Eq. (33) and Eq. (35) as shown in Fig. [5j While 
residual errors of 20%-40%remain, they are largely con- 
tained in a 7r/2 out of phase component rather than in 
an overall amplitude error, and are related to the non- 
linear analogue of the corrections discussed in Appendix 
[B] Furthermore, the modulation of the oscillations comes 
from the power spectrum prefactors of the analytic ex- 
pression causing the inner and outer envelope effect seen 
in Fig. [5] In the signal-to-noise calculation that follows, 
such a modulation mainly cancels out and in fact the use 
of the simple unmodulated C — >■ form of Q jC in Fig. [4] 
suffices out to C ^ 1 for our purposes. 



CMB POWER SPECTRUM 



In |V A| we discuss the phenomenology of a sharp step 
in the inflaton potential on the CMB temperature power 
spectrum and derive scaling relations based on projec- 
tion effects. We then consider WMAP constraints on the 



location and height of the step in S V B 



A. Scaling relations 

The oscillatory features in the curvature power spec- 
trum transfer onto the temperature anisotropy spectrum 
in a manner that reflects projection onto the recombi- 
nation surface, evolution through the acoustically oscil- 
lating plasma, and gravitational lensing after recombina- 
tion. These same effects impact the angular bispectrum 
as well and so it is useful to obtain some physical intu- 
ition in the simpler case of the power spectrum. 

Projection effects damp the amplitude of the fc-space 
oscillations for modes where kr\f ^1. It is simple to 
derive scaling relations for this effect in the flat-sky ap- 
proximation. Ignoring for the moment the acoustic evo- 
lution by taking the Sachs- Wolfe limit, the temperature 
field a(n) = AT/T, in the angular direction n, is given 
by 



1 



a(n) = -p7^(x = Dn) 



1 /■ d^k 



5 J (27r) 



7^ke'''■'^", (64) 



where D is the distance to recombination. The power 
spectrum then becomes 



C,, = 



^ J^PnO^={l/D,k\\)), (65) 



where |j is the direction along the line of sight, orthogonal 
to the plane of the sky. 

Taking the ki]f ^ 1 and d — > limits of the 
change to the power spectrum in Eq. ( 32 ) and assum- 
ing g wconst. so that 



27T A2 



K,0 



p 52 



(66) 
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FIG. 7. Fractional difference of the temperature power spectrum 
(relative to a model with c = 0) of a model with c = 10~^ and d = 
10~^ and different step positions. Gravitational lensing is taken 
into account in these examples. Note that as D/r)f — > 1 (in the 
upper panel), the nodes shift, and the shape changes significantly. 
In black dot-dashed lines we show the expected envelope 
from Eq. ItoI. 



the fractional change in the power spectrum is given by 
AC, 



Cf 



(67) 



In Fig. [6] we show a comparison between the full nu- 
merical result and these expectations for a model with 
Vq — TO^0^/2, and parameters defined in Table [ij With 
this model, rjf = l.M Gpc and D = 14.18 Gpc. We take 
for illustrative purposes c — 10~^ and d — 10"^. 

There are two notable differences between the full re- 
sult vs. the flat sky Sachs- Wolfe scaling. The first is that 
the oscillations are modulated by the acoustic transfer. 
This reflects the fact that at nodes in the acoustic oscil- 
lations between the acoustic peaks there is no transfer 
of power to local temperature fluctuations. The second 
notable effect is a stronger damping starting at £ ^ lO'^. 
This is due to gravitational lensing as demonstrated in 
Fig. [Hj Without lensing the envelope follows the £'^^'^ 
scaling of Eq. ITO] ) as expected. For this reason, when 
considering the bispectrum where lensing effects are more 
difficult to calculate, we will always take an fmax = 2000 
where lensing effects become order unity for the power 
spectrum oscillations. 

It is also interesting to explore the scaling with the fea- 
ture scale rjf. For a model with features at the WMAP 
power spectrum glitches ??/ ~ 1.44 Gpc. For changes by 
a factor of a few around this value, the results scale as ex- 
pected as shown in Fig. [7j However, as the feature scale 
approaches the distance to recombination, projection ef- 
fects take on a very different character. As discussed in 
Appendix [Xj sky curvature prevents a superhorizon fea- 
ture scale from leaving any imprint on the CMB. 

For this reason, when considering the bispectrum 
where we calculate in the flat sky approximation, we re- 
strict ourselves to 77^ < 10 Gpc. 



with the projection factor 



Fix) 



: cos(xz) 



cos(x -I- 7r/4), 



(68) 



where the approximation is in the a; 3> 1 or £ ^ D/irjf 
limit and can be proven by considering that dPP/dx^ = 
(7r/2)Yo(a^)- Rapid oscillations in fc-space are therefore 
suppressed in ^-space by a factor of (£/£/) where £/ = 
D/rjf. It is useful then to scale the numerical results to 
this expectation by defining the scaling factor 



Ar 



\lGvcJ 



-1/2 



(69) 



Note that dividing out by the factor of Ac appropriately 
rescales the oscillations to be order unity at low £ and 
damp as according to the envelope 



2 \D/{lGpc) 



-1/2 



(70) 



B. Constraints 

The WMAP7 temperature and polarization power 
spectra place constraints on the height, width, and loca- 
tion of a step in the inflaton potential. These constraints 
then limit the observability of corresponding features in 
the bispectrum. In order to limit models to reasonable 
cosmologies we also add the following data sets: Bl- 
CEP and QUAD, which include polarization constraints 
dilig, UNION^ supernovae data, the SHOES mea- 
surement of Hq = (74.2 ± 3.6) km/s/Mpc [30] and a big 
bang nucleosynthesis constraint of ilbh^ = 0.022 ± 0.002 
pr , . We include the effect of gravitational lensing on the 
CMB. 

With these data sets and their likelihood functions, we 
perform a MCMC likelihood analysis on the joint step 
and cosmological parameters. We parametrize the initial 
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curvature power spectrum as 



InA^(fc) = liiAs + (ris - l)ln 



Ac f Vf ^ 



O.OSMpc- 



(71) 



3 VlGpc 



V 



W'{kr]f). 



Since we are interested in the Xd ^ oo limit, we take 
a sufficiently small width d in which the damping be- 
havior falls outside of the WMAP range of observation. 
The initial conditions are then described by 4 parameters 
{As,ns, Ac,hirif}, to which we add 4 cosmological pa- 
rameters {flbh'^ ,Qch^ ,9,t} in a flat ACDM context. We 
use Ac as the normalization parameter since we expect 
its errors to be roughly independent of ?]/. 

We take flat priors on each of these 8 parameters. For 
all but ln?7j the prior range is an uninformative > 45cr 
of the posterior in each. We start with a coarse analysis 
in a wide range in 0.4 < rjf/Gpc < 12. Errors on Ac 
throughout the whole range are approximately a{Ac) 
0.03 — 0.06, consistent with the scaling arguments above, 
but for certain discrete ranges of 77/ there is a preference 
for non-zero mean values. These multiple maxima make 
a global MCMC analysis highly inefficient. 

In Fig.jSj we instead show the result of separate MCMC 
chains at fixed values of 77/ away from these special re- 
gions at widely separated r/f = 1.44, and 11 Gpc. The 
results here are consistent with Ac = and yield 1-sided 
95% CL bounds of Ac < 0.05,0.10, respectively. 

Around specific values of rjf a model with Ac ~ 0.1 
is actually a better fit to the data than Ac — 0. In 
Fig. [8] we show an example with rjf — 8.1 Gpc. Pref- 
erence for high frequency oscillations in the power spec- 
trum around the first acoustic peak in the WMAP data 
have previously been noted by ref. |32| and more recently 
by refs. [33t [34] in different contexts. In Fig. |9j we plot 
an example from the 2D chain. Here ?7/ = 8.163 Gpc 
and Ac = 0.11. Note that this model improves the to- 
tal likelihood by 2AlnL = 11.5 for the two extra pa- 
rameters. Moreover, this improvement comes from the 
WMAP likelihood with 2AlnLwMAP = 11-6. Impor- 
tantly, it does not come from the low multipole moments 
nor the glitches at £ = 20 — 40 but mainly oscillations 
with a period of A£ « 5.4 around the first acoustic peak 
100 < i < 300 with contributions continuing out to the 
beam scale. Small variations in rjf around this value 
also produce similar improvements as long as the angu- 
lar scale r]f/D remains fixed via allowed adjustments of 
other cosmological parameters, mainly Q,mh?. 

Other specific values of rjf show comparable likelihood 
improvements, for example at 77/ = 5.135 Gpc and Ac — 
0.09, 2AlnL = 9.9. While this level of improvement for 
2 parameters is notable, it is possible that it represents 
fitting of excess noise at the few percent level in Ci that 
is not well modeled in the likelihood function. 

If these improvements are in fact signal and not noise, 
then there are sharp predictions that can be verified with 
future data. We shall see in the next section that the 
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FIG. 8. Posterior probability distribution of the oscillation am- 
plitude Ac in the angular power spectrum for representative value 
of r)f = 1.44, 8.1 and 11 Gpc. The constraint weakens as r]f ap- 
proaches the horizon for typical cases. For r)f ^ S.lGpc and other 
specific values, the peak of the posterior is shifted to Aq ~ 0.1 
with comparable distribution width. 



bispectrum is observably large so long as the oscillations 
persist undamped to £ ^ 500. Furthermore there must be 
matching features in the polarization power spectra (see 
Fig.|9|. The oscillations in the polarization power spectra 
are comparable but larger at peak as they are less affected 
by projection and modulated by an orthogonal acoustic 
transfer. Detection of this matching signal in polarization 
would make a convincing case for a primordial origin of 
the improvement. 



VI. CMB BISPECTRUM 

We can now use our analytic results to estimate the 
signal-to-noise ratio (SNR) in the bispectrum. We review 
the SNR calculation for a cosmic variance limited data set 
VI A| and approximate methods for its estimation in 



VI B From this approximate method, we derive scaling 



relations for its dependence on the amplitude, width and 
location of the step feature in § |VI C and calibrate them 
against numerical calculations in § VID 



A. Cosmic variance 

The temperature or angular bispectrum is defined as 
the three-point function of the spherical harmonic coeffi- 
cients aim of the temperature anisotropy 
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FIG. 9. Fractional angular power spectrum differences between 
the maximum likelihood sharp-step model {Ac = 0.11, ?;/ = 8.163 
Gpc) and pure power law model [Aq = 0). The WMAP data pre- 
fer a few percent oscillation in the temperature power spectrum 
(top panel) around the first acoustic peak by 2AlnL = 11.6. This 
model predicts matching B-polarization power spectrum oscilla- 
tions (bottom panel) modulated by an orthogonal acoustic phase. 



The cosmic variance of the Gaussian part of the field puts 
an irreducible limit on the SNR of 



C e,Cp^Ci>„de, f„f„ 



(73) 



where 



accounts for permuted contractions of repeated ^'s and 
the angular power spectrum is defined by 



(75) 



B. Approximations 

Here we are interested only in an order of magnitude 
estimate for the SNR through a crude computation of the 
angular bispectrum from the curvature bispectrum. We 
therefore take the flat-sky approach and the Sachs- Wolfe 
limit for the temperature anisotropy. 

In the flat sky approximation, the angular bispectrum 
is defined by the three-point function of the Fourier mo- 
ments of the temperature field given by a(l) 

{a{h)a{h)a{h)) = (2^)H{\i +h + h)B(^,^,e,.e,y (76) 

For £i, £2, ^3 ^ Ij it is related to the all-sky bispectra as 



B, 



(2£i + 1)(2£2 + 1)( 24 + 1) / 4 £3 

I 



(77) 



Under the flat Sachs- Wolfe approximation of Eq. ([64| , 
the bispectrum becomes 



B, 



where 



2 dfciii 
la 



2tt 



-Bn{ki,k2,k'i), 
(78) 



ki = (li/i?,fci||), k2 = (l2/i?,fc2||), k3 = -ki-k2. 

(79) 

It is useful to note that with the correspondence be- 
tween all-sky and fiat-sky expressions |35| 
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(80) 



Scaling relations for the SNR can thus be derived from 
scaling relationships for the flat sky spectra combined 
with a counting of triangles in the available phase space. 



C. Scaling arguments 



We thus require an efficient means of predicting the an- 
gular bispectrum given the curvature bispectrum. 

The bispectrum arising from the approximations in 
Eq. (48) and Eq. (29) looks especially formidable to 



Before we numerically compute Eq. (73), we can esti- 



project onto the angular sky due to the inseparable 
damping function T>{x). In Appendix [C| we demonstrate 
that, to a good approximation, this bispectrum can be 
cast into an approximately separable form, which will 
make the full projection onto the angular sky a much 
more tractable problem. Still, this separable form re- 
quires computational intensive operations making an ex- 
ploration of the whole parameter space difficult. We leave 
evaluation of this to future work. 



mate how strongly we expect the SNR to scale as we add 
triangles. 



1. Local bispectrum 

As a proof of technique, let us first examine the case 
of a local non-Gaussianity where the reduced bispectrum 
is a constant or 

Bn{ki,k2,k3) = ^/nl [Pn{ki)PTz{k2) + perm.] . (81) 
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The signal is dominated by squeezed configurations is ^ 
£l and the above arguments imply that 



nlocal 



{1^1 D) Al {Is/D) 



(82) 



Consequently, the SNR for triangles with long and short 
sides between imin and £niax for a scale invariant power 
spectrum will go as 



In 



(83) 



where we have used Eq. ( 66 1 to take Ci (x I in the 



same flat-sky, Sachs- Wolfe limit. 

The integrals can be thought of as counting triangles: 
there are £| ways of choosing one of the long sides of the 
triangles and likewise ^| ways of choosing the short side. 
The final long side is determined by requiring the triangle 
close and so the total number of triangles becomes i'j^i'g- 

The naively infrared divergent integral is regulated by 
the lowest multipole available, which in our cosmic vari- 
ance limited calculation is imin = 2. This leaves an over- 
all scaling of ^^^^ i in agreement with the cosmic variance 
limit of a well known result in the literature that the er- 



ror on /nl drops as N 



-1/2 
pix 



where TVpix = /sky^max is the 



number of observed pixels |36l I37j . 



The efi^ect of the damping envelope ^{x) is to impose 
a limit on the maximum scale that can contribute to 
the SNR, ^max- Neglecting the prefactors and taking 
£f — D/rjj, we expect the projection of the bispectrum 
in Eq. ( 84 ) onto the flat-sky to scale as 



1 
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(85) 



where the [KrjfY supplies two factors of iijlj) and the 
two integrals over fcy suppress the result by {ilif)~^ due 
to the oscillatory integrands as in the power spectrum. 
Finally, the integrals over l\ and li provide a factor of 
and yield 
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The ^max factor can be thought of as a counting of equi- 
lateral type triangles: in three dimensions, there are ^^ax 
ways to choose the first side, another ^^^x ways to choose 
second side while the third side is determined by requir- 
ing the triangle close, yielding ~ f^^x triangles. 



Flat type 



2. Feature bispectrum 



Flat triangles scale equivalently to equilateral type tri- 
angles. 



We can now look at our feature bispectrum. There 
are now three cases we need to consider, the SNR being 
dominated by contributions from: 

1. equilateral type shapes: ki ^ k2 ^ k^] 

2. flat type shapes: /ci ~ ^2 ^ 2/03; and 

3. squeezed type shapes, as in the local case k^ <C 
ki ^ k2. 

We consider these separately below, estimating how their 
contribution to the SNR scales with the number of 
modes. 



Equilateral type 



Consider the leading contribution to the bispectrum of 



Eq. (48) assuming a scale invariant spectrum, and taking 



/q^— Afi{k). In the UV or large-fc limit, the bispectrum 
is dominated by the term quadratic in the perimeter of 
the momentum triangle 
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however, due to the restriction that all modes be co- 
linear, this reduces the number of triangles available. 
Flat triangles correspond to the restriction ^2 — ^sl < 
AL, where L = £i+£2+i3 and AL is the tolerance. With 
this restriction, we expect the number of triangles to scale 
^s ^max V ^-^/^max- The curvcd-sky imposes a minimum 
AL = 1, and thus we expect the number of triangles to 
scale as -^^ax and the contribution of flat triangles to the 
SNR scales as 



/>3.5 
^max 



flat 



«-2 /)-6 
/ ^max 

■tmax 



/-2/)3.5 
^/ ^max- 



Thus flat triangles contribute less to the total SNR than 
equilateral triangles for ^max ^ £f- 



Squeezed type 

Considering now squeezed triangles, in the limit ka, 
^L, and taking the leading order term 



Bn{kL,kL,ks)~ - -i2n^ ^ 



Xd 
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(89) 
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we estimate that this bispectrum contributes to the signal 
in squeezed angular space configurations is ^ as, 



R /)-5/2fl-3/2 



(90) 



Compared to the result in Eq. (85 1, for a fixed (.q, the sig- 



nal in squeezed triangles falls slower in the UV relative to 
the signal in equilateral triangles. However, considering 
the SNR, 



N 



we obtain 



g3 
max 



(91) 



(92) 



which implies that the contribution of squeezed triangles 
compared to the equilateral type goes to zero as ^max — >■ 
oo. Interestingly, while squeezed triangles dominate the 
signal at high £, they do not dominate the signal to noise. 
The reason for this is two-fold. Squeezed triangles in 
themselves suffer from higher cosmic variance which, for 
a given triangle, eliminates the scaling advantage. In 
addition, there are many less triangles which contribute 
to the squeezed limit. 




FIG. 10. The square of the SNR per unit bispectrum ampli- 
tude (Ab) for id = 8514 for all triangles, and for flat config- 
urations only. We also plot the scaling relations from j ]VI C| 
In the lower panel, we demonstrate that changing rjf while 
fixing Ab and Id leaves the SNR fixed at high £. 



D. Scaling and numerical results 

To summarize the scaling relations of the previous sec- 
tion, for small values of C, 



N 



12 ^max _ q2 max 



JJ2 ' 



(93) 



as c? ^ and for fmax ^ With this behavior in mind, 
we define the angular bispectrum amplitude. 



Ah=C 



IGpc 
Vf 



6c 



IGpc 



3/2 



3c \lGpc 



(94) 



which takes into account the dependence of the amplitude 
on the step scale, rjf. Compared with Ac the analogous 
amplitude for the angular power spectrum, the bispec- 
trum amplitude increases as rj^^. Recalling that the ob- 
servational errors on Ac are only weakly dependent on 
77/, the SNR is maximized by placing the feature at the 
largest observable scales, near the current horizon. 

With a finite d there is an effective maximum multipole 
beyond which the SNR saturates. This scale depends on 
77/ as 



2D 2D V2eo + 6c 

a^rf = 1 ■ 

j]f rjf TTcL 



(95) 



In what follows we simply scale out the multipole space 
bispectrum amplitude Ab, and present the square of the 
SNR per unit as a function of the multipole space 
damping scale, £d or ^max- These redefinitions render our 
results independent of the choice of the model parameters 
{(/)/, c, d}. Of course to make connection with a particular 
model, one must rescale by the appropriate factors oi Ab 
and choose the correct damping scale. 

We verify the SNR scaling behaviors and obtain the 
proportionality coefficient in the scaling relations by nu- 
merically evaluating Eq. (73). In the upper panel of 
Fig. [To] we show the contribution to the SNR of all trian- 
gles and the contribution from triangles that are exactly 
flat. The close agreement of i"^ and £^'^ curves with the 
numerical results verifies our scaling relations, confirm- 
ing that the SNR is dominated by triangles whose sides 
are all of comparable length. For this test, the feature to 
occurs at a scale rjf = 1.44 Gpc. 



The lower panel in Fig. 10 demonstrates that our re- 
sults are independent of the specific choice of 77/; plotted 
is the square of the SNR per unit bispectrum amplitude 
under a shift of the feature from 1.44 Gpc (black, solid 
curve) to 2.88 Gpc (red, dashed curve), while simulta- 
neously reducing the width of the step as to hold the 
damping scale, £d, constant. 

In Fig. 11 we show the dependence of the SNR on 
the angular damping scale id- Notice the SNR saturates 
once the maximum multipole ^max exceeds the damp- 
ing scale. In fact, to excellent approximation, the value 
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at saturation can be approximated from the d = re- 
sults by setting ij^ax = l-06£d- However, for values of d 
where id > 2000 secondary anisotropy contributions such 
as gravitational lensing and the Sunyaev-Zeldovich effect 
become important. Furthermore the Planck sateUite will 
be limited by its instrumental beam to lower multipoles. 

Thus our rough approximation of the maximal signal 
to noise accessible to the CMB temperature bispectrum 
can be expressed in terms of our scaling parameters and 
the amplitude of the features in the temperature power 
spectrum as 

where £raax — min(1.06£rf, 2000). Combined with the up- 
per limit of Ac ^ 0.05 for a typical value of 77/, (see 
Fig. [s]), one infers that sharp steps with 77/ > 2 Gpc have 
potentially observable features in the bispectrum while 
still being consistent with current power spectrum re- 
sults. For the maximum likelihood model, where Ac = 
0.11 and 7]f = 8.163 Gpc, the maximal (S/N)^ sa 790. 
In fact even taking ^^ax — 800, the {S/N)"^ ^ 20 and so 
there may even be information about this model in the 
WMAP data if an optimal analysis can be performed. 

Our scaling relation makes it seem that we can in- 
crease the SNR without bound by moving the feature 
to larger and larger scales. However, just as was the 
case in the power spectrum, our flat-sky approximation 
overestimates the bispectrum as 77/ approaches the cur- 
rent horizon due to projection effects on the curved sky. 
Adopting the same maximal scale of rjf < lOGpc as was 
found for the power spectrum, the SNR can be at most 
comparable to that of the maximum likelihood model. 

Should a horizon scale feature be detected in the data 
then one should also consider the sample (co)variance 
of the finite number of such features that can fit in our 
horizon volume when measuring the parameters associ- 
ated with the feature. For upper limits, the Gaussian 
noise variance suffices. 



VII. DISCUSSION 

We have investigated the power spectrum and bispec- 
trum of curvature fluctuations generated by a sharp-step 
feature in the inflationary potential. In the limit where 
the change in potential energy due to the step is small 
compared to the kinetic energy of the inflaton, we have 
demonstrated that the background system can be solved 
perturbatively in this small parameter. 

Using this solution for the background evolution of the 
inflaton, we obtained closed form, approximate analytic 
solutions for both the power spectrum and bispectrum of 
curvature fluctuations. The analytic solutions generically 
take the form of an oscillating window function times a 
damping function. The form of the window function is 
set by the behavior of the step in the inflnitely sharp limit 
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FIG. 11. The square of the SNR per unit bispectrum ampli- 
tude (Ab) for various damping scales Id- The shaded regions 
indicate the prediction of our scaling formula in Eq. ( |96| l and 
match the asymptotic high £ ^ Id saturation point for each 
case. 



with its characteristic frequency determined by the step 
position in conformal time, ijf. The damping function 
on the other hand determines the maximal wavenumber 
out to which the oscillations persist, which scales with 
the inverse of the finite width d, with a form factor that 
depends on the shape of the step. In this work we have 
explicitly calculated the damping function for a step with 
a hyperbolic tangent shape. 

Our analytic power spectrum and bispectrum solutions 
pass several important consistency checks. They van- 
ish for modes that are superhorizon size when the field 
crosses the step, and thus do not introduce spurious su- 
perhorizon effects. They satisfy the consistency relation 
which relates the squeezed limit of the bispectrum to the 
tilt of the power spectrum, or scalar spectral index. Ad- 
ditionally, leading order corrections to this relationship 
scale only as the square of the squeezed side. Numeri- 
cally, for very small perturbations to the kinetic energy, 
our solutions are accurate to approximately 10%. This 
difference is largely due to slow roll corrections to the 
mode functions, which introduce additional phase shifts 
and thus alter the form of the non-Gaussianity. Some- 
what surprisingly, rather than one percent level correc- 
tions as one might naively guess, these type of corrections 
generically lead to 10% corrections to the bispectrum. 

As the step height increases, and the perturbation to 
the kinetic energy becomes large, the analytic solutions 
begin to break down due to two effects. First, the per- 
turbation to the kinetic energy becomes large and infla- 
tion ends at an appreciably different time for a potential 
with and without the step. This difference corresponds 
to a change in the matching of inflationary fleld scales to 
physical scales and hence a phase shift between the an- 
alytic solution and its numerical counterpart. We make 
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use of the approximate time translation invariance of the 
inflationary background to absorb this phase shift into 
a redefinition of the conformal time of the unperturbed 
background. The second effect is that as the kinetic en- 
ergy perturbation grows to order unity, our analytic so- 
lutions attain a larger amplitude while damping away at 
a faster rate compared with the numerical results. In our 
analytic solution these properties are controlled by a sin- 
gle parameter, e^f, which effectively measures the ratio 
of kinetic to potential energy of the inflaton. At leading 
order, our analytic solution takes into account only the 
kinetic energy of the inflaton on the background without 
the step. However, we have shown that the discrepancy 
between the numerical and analytic results can be sub- 
stantially improved by including the first order kinetic 
energy perturbation in e^- 

As the step width becomes infinitely sharp, the high 
wavenumber limit of the spectra exhibits some interest- 
ing behavior. We have found that the perturbation to the 
curvature power spectrum approaches a constant ampli- 
tude oscillation, while the reduced bispectrum, or effec- 
tive non-linearity, diverges quadratically in the perime- 
ter of the momentum space triangle. The reduced bis- 
pectrum or non-linearity parameter /nl can thus reach 
10'* or more. While one might naively infer a strongly 
divergent non-Gaussianity on small spatial or angular 
scales, the real space analogues of the power spectrum 
and bispectrum, the two- and three-point correlation 
functions remain tiny for most configurations, in partic- 
ular all those involving only small separations. Oscilla- 
tions largely prevent the various triangle configurations 
adding up to anything significant: For the most part they 
simply average to zero. However, for certain configura- 
tions, where the separation between points is compara- 
ble to the physical scale of the feature, resonance in the 
integral leads to logarithmically large values of the cor- 
relation. High wavenumber oscillations in the spectra 
therefore are the result of Fourier transforming a sharp 
feature in the correlation functions on large scales. 

The real space picture of the bispectrum suggests that 
the fc-space triangles are not all independent as they must 
sum to yield a sharp feature at a fixed, large scale. This 
covariance of triangles must be accounted for if such a 
signal in the bispectrum is detected in the future. In prin- 
ciple, one could use our methods to evaluate higher Ap- 
point functions which quantify this covariance and that 
of the power spectrum. Furthermore, they would pro- 
vide corroborating evidence for a step and are expected 
to grow even more strongly with the momentum perime- 
ter than the bispectrum considered here. 

On the other hand, to assess the detectability of the 
bispectrum, it is sufficient to compute the SNR where the 
noise is attributed to Gaussian random fiuctuations only. 
We utilize our analytic bispectrum solutions to assess de- 
tectability as a function of the step location, height and 
width. For simplicity, we neglect acoustic transfer and 
the sky curvature, working in the Sachs- Wolfe and flat- 
sky limits. This over-simplification suffices for our order 



of magnitude estimate. Under this approximation, the 
signal-to-noise ratio is dominated by equilateral trian- 
gles and scales with the fourth power of the maximum 
multipole due to resolution or finite width d. These scal- 
ing relations allow us to characterize observability as a 
function of the step parameters. 

The SNR is maximized by placing a sharp feature close 
to the horizon scale. Since the effective non-linearity 
grows quadratically with the ratio of the feature scale 
to the observed wavelength, placing the feature at larger 
scales means it can grow to a larger amplitude within 
the window probed by the CMB. The maximum scale at 
which one can place the feature before it becomes unob- 
servable is the horizon scale for both the power spectrum 
and bispectrum. In real space, the high-fc behavior of 
either translates to sharp correlation features at separa- 
tions comparable to the physical scale of the step. Be- 
cause of causality, an observer confined to a single posi- 
tion cannot measure a correlation across a distance larger 
than the horizon. Since in the d — limit oscillations in 
the power spectrum persist even for superhorizon scale 
features, the momentum space results naively seem prob- 
lematic. However, as the horizon scale is approached, the 
frequency of the oscillations in the spectra approach the 
fundamental spacing of the spherical harmonics, A£ = 1. 
Thus, as this scale is exceeded, CMB spectra are ren- 
dered insensitive to the presence of the feature, except 
for an unobservable shift in the monopole. 

While one might guess that the angular power spec- 
trum of the CMB would place severe restrictions on large, 
sharp steps due to their oscillations persisting into the 
strongly constrained £ ^ 200 — 400 region, in fact the 
constraints are considerably weaker. There is an addi- 
tional source of damping when the oscillating features in 
the momentum space power spectrum are projected onto 
the sky. Additionally, effects such as gravitational lensing 
tend to hide the oscillations on very small scales. Curi- 
ously, the WMAP power spectrum data are in fact better 
fit by such a feature than a pure power law spectrum with 
2AlnL « 12 for the 2 extra parameters. While this im- 
provement may be due to fitting excess noise in WMAP, 
it is formally more significant than a similar fit of the 
well-known glitch at ^ = 20 — 40 with 3 extra parameters 

If this improvement is due to a slow roll violating fea- 
ture such as a step, there are testable consequences in 
both the polarization power spectrum and the bispec- 
trum. The polarization power spectrum should carry 
matching oscillations whose amplitude is less affected by 
projection. The bispectrum would be detectable as long 
as the oscillations continue undamped to at least i ^ 500. 
In this case, the bispectrum may be able to confirm the 
primordial origin of such a signal. 
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Appendix A: Real Space Correlation 



Given that the step feature is localized in time, it is 
interesting to consider the correlation functions, the real 
space analogs of the power spectrum and bispectrum. 
While somewhat orthogonal to the main thrust of the 
paper, we will see that the real space correlation func- 
tions give insight into causality, sample vs. noise vari- 
ance, consistency relation and the perturbative validity 
of our expansion of the action. In particular, we shall 
see that the k ^ oo behavior of the power spectrum and 
bispectrum correspond to sharp features in the real space 
correlation functions. 



In SAl 



This Appendix is organized as follows 
demonstrate that, in the limit in which d = 0, the per- 
turbation to the curvature power spectrum leads to a 
sharp feature in real space, the modes as fc — >■ oo sum- 
ming to give a logarithmically divergent derivative at a 
single point associated with the physical scale at which 
the feature occurs. In |A 2[ we show that, despite attain- 
ing a large non-linearity, /nhik), the real space analogue 
of the bispectrum remains small everywhere, except in 
the vicinity special points, associated with the physical 
scale of the feature, where it becomes large, diverging as 
~ ln{d). This behavior is shown to be due to the fact 
that the momentum space triangles only coherently add 
for very special configurations of the spatial points which 
lead to resonances in the Fourier integrals. In |A 3[ we 
show that at least some of this behavior can be antic- 
ipated from the consistency relation, which relates the 
slope of the two-point correlation function to a certain 
limit of the three-point correlation function. 



The two-point correlation function is the Fourier trans- 
form of the power spectrum 

6K(ri,r2) = (7^(rl)7^(r2)) 
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k kr 

where r = |ri — r2|. 

With our analytic calculation for the power spectrum 
of the step feature in |III| we showed that in the limit in 
which the step width d — )■ 0, the power spectrum correc- 
tion attains a constant amplitude oscillation which per- 
sists to fc — >■ oo. Corresponding to this high frequency 
behavior, there must be sharp features in the correlation 
function. 

The contribution of the step feature to f (r) can be 
calculated from Eq. ( 32 ) as 



dk sin(fcr) 
k kr 
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When d = 0, Eq. ( A2 1 can be evaluated analytically, and 
one obtains, 
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(A4) 



We plot this function in Fig. [12] In this case, the sum 
over all of the modes out to infinity results in a sharp fea- 
ture in the slope at r = 2r// while the function remains 
finite, and order 0{C) everywhere. The fact that r = 2rif 
is a special point can in fact be read directly off of the 
integrand. Since W oc cos{2kT]f) at high k and the expo- 
nentials in the transform give sin(/cr), contributions from 
different fc-modes oscillate away except for the stationary 
phase point r = 2rif. We will use this type of reasoning to 
deduce the behavior of the three-point function in §A 2 a| 

The slope of the correlation function £,2n{r) formally 
diverges for d = at the point r = 2rjf. More generally 
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and this integral is logarithmically divergent with d at 
r = 2r]f. We will relate this divergence to the bispectrum 
behavior in fjXSl 



Examination of Figs. [T] and 12 also raises an interest- 
ing question about causality. As we have demonstrated, 
a sharp feature in the inflationary potential gives rise 
to oscillating features in the power spectrum of curva- 
ture fluctuations and as the step becomes sharper and 
sharper, the oscillating features ring to larger and larger 
wavenumbers, k, before decaying. One might then be 
tempted to conclude that an observer who only had ac- 
cess to a region of space 77 < 77/, could measure superhori- 
zon perturbations, by observing oscillations in his power 
spectra on small scales. 

We have, however, also shown that, in real space, these 
modes simply sum up to a sharp feature at the physi- 
cal scale corresponding to the size of the physical hori- 
zon when the inflaton crosses the feature. One therefore 
comes to the conclusion that an observer who has access 
to only a region of space 77 <C ?// cannot measure the fea- 
ture. The resolution to this apparent discrepancy is that 
the fc-space features are spaced by less than l/rj and so 
cannot be sampled in the finite volume. In the CMB this 
is imposed by the finite volume within the distance to 
recombination. In large scale structure, this is similarly 
imposed by the finite survey volume. 

More specifically, in the case of the CMB the projec- 
tion of these oscillations onto the angular sky imposes a 
very physical cutoff on the frequency of the oscillations, 
corresponding to the fundamental spacing of the multi- 
poles, A£ = 1. Oscillations in the power spectrum with 
a frequency spacing smaller than this fundamental are 
rendered unobservable when projected onto the spheri- 
cal sky. In the Sachs- Wolfe approximation. 
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In the example at hand, a,t krjf ^ 1 



A2,(fc) « A|j_o(fc) [l-Ccos(2fcr;y)] 



(A8) 



Since the Bessel function oscillates only at a frequency 
kD, if 277/ ^ D then the oscillatory piece cancels out of 
the integral. This cutoff prevents an observer from ob- 
taining information about sharp features on scales larger 
than their horizon. However, we note that this cutoff 
is only exactly enforced by the full projection onto the 
spherical sky. In particular, in the periodic flat-sky ap- 
proximation despite enforcing a fundamental spacing of 
A^ = 1, features with oscillations at A^ < 1 can im- 
print beat frequencies on the re mai ning modes making 
SCi/Ci. (X (irjf/D)-'^^^ as in Eq. Fig. [t] verifies that 
this scaling breaks in the full sky calculation ioT rjf > D 
leaving a suppressed contribution to the anisotropy. 



0.2 



0.1 



<] 

S3 



-0.1 



-0.2 - 



-0.3 - 




FIG. 12. Real space two-point correlation function changes 
due to the addition of a step feature (d — ^ 0). At r/ry/ = 2, 
the slope of the correlation function diverges as the dashed 
line indicates. 



2. Three-point correlation 

Similarly, the three-point correlation function is re- 
lated to the bispectrum by a Fourier transform of its 
arguments 



C3K(ri,r2,r3) = (7^(rl)7^(r2)7^(r3)) 
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For a general configuration and the step potential bis- 
pectrum, this integral is difficult to compute. As was 
the case with two-point function, we can gain insight on 
its behavior by first studying resonances in the integrand 
where the phase is stationary. We can then confirm this 
behavior by explicitly evaluating simple configurations. 



a. Resonances 

The bispectrum of the step feature as its width rf — > 
and the wavenumbers ki — )• 00 have oscillatory behavior 
given by Eq. (l56|) 



B-iziki, k2, ks) cx cos[(A:i + k2 + k3)r]f] 
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This implies that in the real space correlation function 
most of the contributions will integrate away except for 
special resonant points where the frequencies in the trans- 
form match those of the bispectrum. 

Without loss of generality, we can set ra = and inte- 



grate out ks in Eq. ( A9 1 



67^(ri,r2,0) = ny ^e^'^'-'i37^(fci,A;2,fc3), (All) 
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FIG. 13. Relationship between bispectrum triangle shapes 
and resonances in the three-point real space correlation func- 
tion. As d — >■ each set of self-similar bispectrum triangles de- 
termined by (ki, k2, k3)s, where s is some scale factor, causes 
a resonance at (ri,r2,r3) and its mirror image. As the over- 
all size of the triangle s — >■ oo this resonance becomes sharper 
and sharper in real space but remain associated with sepa- 
rations 0{rjf). Flat triangles are responsible for the largest 
scale resonance for ri — ra = r2 — ra — 2rjf. 



where ^3 is now defined by the cosine law. The inte- 
grals are thus sums of bispectrum triangles defined by 
the 6-dimensional vector K = (ki,k2). In this 6D vector 
notation, the phase function 
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= K'^R± [yK^tiK + VKTlaK + VKTUKJ 7?/, 
where R = (ri,r2) and we have defined 
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Now, we can look for points where this is stationary with 
respect to triangle configuration. Setting d(j)/dK. ~ 0, we 
obtain 
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which can be written, 

= T»7/(kj - kg). 
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where hats denote unit vectors. Note that the real space 
resonances occur for triangles that are coplanar with the 
fc-space triangles. At these solutions, the phase function 
reaches 
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T(fci + fc2) ± (ki + k2) • k3 ± (fci + fc2) ± A:3 
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= 



and so satisfies the constancy condition. 

For each fc-space triangle geometry there is a resonance 
in real space triangles that is fixed by the fc-space geome- 
try (see Fig. 
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The orientations of relative to kg are 
determined by the half angles of the respective fc-space 
triangle angles. The physical length varies from a max- 
imum of 27yj, if ki II — ks; through y/2rjf, if they are or- 
thogonal; to 0, if they are parallel k; || ka. Note, however, 
that if ri — then r2 2r]f so that resonant triangles 
always have characteristic size rjf. The maximum size of 
the correlation region of 2rif is achieved for fiat triangles. 
Finally, there is no resonance at ri = r2 = rg = since 
3 coparallel /c-vectors cannot close. 

The allowed resonances all correspond to points in- 
scribed on a sphere of radius rjf where rotation in or 
around the kg direction corresponds to a rotation of the 
real space resonance. 

We can also show that these points are extrema rather 
than unstable saddle points by considering the Hessian 
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A similar but more tedious exercise shows that this ma- 
trix is positive semi-definite. Namely for any vector 
Q = (qi,q2) defining a direction in the 6-dimensional 
space of bispectrum triangles with a closure relation 

qs = -(qi + q2) 



Q^HQ 



3 

^ h 



(qi • kj 



> 0. 



(AI8) 



The eigenvalues vanish along the special directions of 
qi II ki. There are two cases where this happens. The 
first corresponds to fixing the triangle geometry while 
changing its overall size. The second is for the family of 
fiat triangles where the sides are all coparallel but the 
ratio of lengths fci//c2 can vary. For these triangles there 
then two flat directions whereas non-degenerate triangles 
have only one. Equivalently, Eq. (AI6) shows for both of 



these cases, the phase function is strictly fixed. 

Given the behavior of trajectories in the space of tri- 
angles, we can estimate the volume of phase space that 
satisfies the resonance condition, and thus contributes to 
the integrals. Except for these special flat directions in 
the 6-dimensional space , the typical width of the reso- 
nances is Ag ~ i/fcT??/. Putting this together, we can 
estimate the degree of convergence or divergence of the 
three-point correlation function in Eq. (All ). For a typ- 
ical real space configuration where ri 7^ r2, there is only 
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one flat direction ~ ki where the triangle size is mul- 
tiplied by some scale factor s 



C37?,(ri,r2,0) cx r]f / dk. 



k-^v 



(X 



max'?/ 



(A19) 



where /cmax oc from the damping function D. 

On the other hand, for the special configuration where 
ri = r2 — 2rif, there are two flat directions, kg and R = 
^2/^1 which implies the three-point correlation function 
will scale as 



6K(ri,r2,0)cxr;2 J k^dk, ) k'^V 
oc ln(fci„ax'y/), 



(A20) 



where we assume that ki > Tj'J^ . 

The Hessian also identifies another special case where 
there is an approximately flat direction if fc^ ^ TyJ ^. This 
is the case of squeezed triangles where one of the sides 
does not contribute significantly to the phase function. 
In this case k^ <^ k2 ~ ki = k 



k • (ri - r2) ± 2fc77/, 
for which the stationary solutions are 
ri - r2 = ±2?7/k. 



(A21) 



(A22) 



such that there is a resonance whenever the separation 
between any two points is 2rif, regardless of the position 
of the third. In this case the dominant term in the bispec- 



trum goes as 77/fc ^fcg ^ [see Eq. (89)] and the resonance 



contains one fiat direction in k parallel to ri — V2 and 
two orthogonal to of width ^fkjrjj each. Therefore the 
three-point correlation scales as 



k-'^v I d^k^k^^ 



6K(ri,ri+2r?/e,0) CX7?/ J dk (^^^ k'^V J 

cx ln(fcmax'?/)ln(fcmin?7/), (A23) 

where e is any unit vector. Here the correlation depends 
on a cutoff in the infrared k^i^ which in practice comes 
from the survey volume of the data, i.e. the current hori- 
zon for a cosmic variance limited statistic. 

These results imply that, despite the naive diver- 
gence of the reduced bispectrum, the three-point correla- 
tion at most positions actually converges except if points 
are separated by 2r]f where fiat and squeezed triangles 
cause a logarithmic divergence. 

The real space picture of the non-Gaussianity of a 
sharp step raises an interesting question as to the mean- 
ing of the noise term in our signal-to-noise calculation. 
In real space, the increasing (S'/iV)^ at high k is as- 
sociated with modes that superimpose to give a sharp 
correlation feature on large scales. Thus each triangle 



cannot be considered an independent probe of the non- 
Gaussianity. By including only the Gaussian contribu- 
tion to the noise, we neglect this sample covariance of 
triangles. Our SNR therefore only quantifies detectabil- 
ity of the non-Gaussianity rather than the measurability 
of parameters associated with the step potential. If in 
the future, non-Gaussianity of this type is detected, then 
the covariance must be considered by computing the 6pt 
functions. 



b. Configurations 

We can test these properties of the three-point function 
by explicitly evaluating it for a few simple configurations. 
The expression in Eq. ( A9 ) involve six integrals over a 



highly oscillatory function, making it difficult to evaluate 
numerically. However, in the limit in which two of the 
points are coincident, the expression can be simplified to 



2 / kidki f fc2*2 



67^(r,0,0) =(2^)^ / ^ / = / k,dk, 
(27r)3 J (27r)3 J^k.-k^l 



X Bn{ki,k2,k'i) 



sin(/cir) sin(A;2'') 
kir k2r 



(A24) 



In the limit d = 0, using the result for the bispectrum in 
\ IV B the integral over k^ in Eq. ( A24 ) can be performed 



analytically (the result is messy and not particularly en- 
lightening, and we omit it here) leaving a more tractable 
problem of evaluating only two integrals. In order to nu- 
merically evaluate the resulting expression, and render 
the result at r/r]f = 2 finite, we need to impose a cutoff 
in the integral at high wavenumber. However, the inte- 
gral needs to be cut off in a way that reflects the fact that 
triangles of different shapes experience different damping 
behavior. For non-zero values of d, we note that as long 
as d is small, this result can be utilized to integrate the 
bispectrum at flnite d, by parts. We then need only keep 
the boundary term, since the remaining term will be sup- 
pressed by a factor of d, and will thus be negligible in the 
limit d — > 0. Moreover, since this explicitly retains the 
form of the damping function, we automatically apply 
the correct damping to each shape. 

As well as divergences at the point r/rjf — 2 arising 
from the ultra-violet parts of the phase space, there are 
also the usual divergences associated with the infrared 
parts of the scale invariant power spectra. In reality, 
these divergences are regulated by the finite size of the 
longest observable wavelength. To perform the numerical 
evaluation, we simply apply a hard cutoff on the mini- 
mum momenta, fcmin- 

To perform the integration, we write the bispectrum 
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1 2 3 4 5 



0.001 0.01 



FIG. 14. Contribution to the real space three-point correla- 
tion function from the three terms in Eq. ( |A25 1. The thick 
dashed curves were computed for damping scale Xd = Xdi = 
200 with an infrared cutoff at fcmin = O.l/rjf, while the corre- 
sponding thin solid curves were calculated with Xd2 = 1000. 



FIG. 15. Behavior of the various contributions to the real 



space three-point correlation function ^snii", 0, 0) in Eq. ( A25 1 
at the point r = 2i]f as a function ultraviolet cutoff Xd- The 
dependence on the infrared cutoff is also demonstrated. Here 
fcmin.i = 0.1/»7/, while fcniin,2 = 0.001/r;/. 



3. Real space consistency 



as 



Bn{k^,k2,k^) « 

1 2 3 

- In{K)kik2k3 + l2{K){kl + kl + ki) 
il2{K)^hiK))Y,k^k,]. (A25) 



In Fig. [14] we show the contribution to the real space 
three-point correlation function due to each of the terms 
in Eq. ( A25 1 at two different values of the step width, 



or damping scale Xd- In particular, note that the skew- 
ness of the field is finite, and small while ^37^(r, 0, 0) van- 
ishes quickly for r/r/f > 2. The three-point function is 
also largely insensitive to the value of the damping scale, 
Xd, which regulates the UV, except at the special point, 
2rif, where it diverges. 



15 we demonstrate 



r — 'irjf, wnere it diverges. In Fig. 
that this divergence is no more pathological than ln((i). 
Fig. [15] al so de monstrates the scaling behavior of Eqs. 
( A20 ) and ( A23 ) , note that the contribution to the three- 
point function due to the most dominant UV term, /q, 
is insensitive to the infrared cutoff, fcmin, while the term 
which dominates the squeezed contribution, I2 — Ii, is 
logarithmically sensitive to the infrared cutoff. The re- 
maining term, in Eq. (A25) is not UV divergent and we 
omit it from the plot. 

All of these properties and scalings are consistent with 
our resonance estimates from 9A 2 a 



To shed light on why the two-point correlation reaches 
infinite slope at r = 277^ while the three-point correlation 
logarithmically diverges there, we can extend the familiar 
fc-space consistency relation between the power spectrum 
and squeezed bispectrum to real space correlation func- 
tions. 

In fact the fc-space consistency argument itself comes 
from a real space derivation. We follow the analysis of 
|25) here but make more explicit what is meant by the 
background. To this end, we split the field up into low 
pass and high pass filtered pieces 



7^(x) ^n + n 



(A26) 



d^k 



k< 



k, (2vr)^ 



k>k. 



d^k 



That is, we split the contributions of the field at each 
point into contributions from modes shorter than some 
reference mode fc*, and long wavelength modes assumed 
to be much longer than this scale. Next we consider the 
two-point function of the high pass field in the fixed back- 
ground. Since the background fluctuation is expected to 
be small, we can functionally expand 



(7^(xl)7^(x2)) =(7^(xl)7^(x2)) 



K=0 

(^(xi)7e(x2)) 



(A27) 



where the background field is slowly varying and evalu- 
ated in the vicinity of the points, which for definiteness 
we take to be x+ = (xi -I- X2)/2. In the presence of the 
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background mode, the metric is 



(A28) 



Then, absorbing the effect of the background mode into a 

^ g7l(x) 

in the separation between points 



change in variable x ^ x' — e^^^'x, we induce a change 



(A29) 



Defining the two-point correlation function of the high 
and low pass fields 



6tc(Axi2) = {nixi)n{x2)) 
dk 



sin fcAcci2 



k>k. 



k '^^ ' fcAxi2 



f27j(Ax+3) = (^(X+)^(X3)) 



dk . n , sin fcAx+s 

— A|,(fc — ^ 

fc<fc. k ^ fcAx+3 ^ 



where Axi 



Ix," — Xil, we obtain 



(7t(xi)7t(x2)) = 6K(Aa^i2) + ^(X4 



27?, 



din Aa;i2 

which then correlates with a third point as 



(A30) 



(A31) 



(A32) 



(7^(xl)7^(x2)^(x3)) - 67j(Aa:+3) 



din Aa^i 



This defines the squeezed contribution to the analogous 
three-point function. 

Note that we can Fourier transform this relation to 
obtain the bispectrum. Given that through integration 
by parts 
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d^k 



dlnr A>fc. (27r)3 
we obtain immediately. 



1 d(fc3p^) 



k^ d\nk 



(A33) 



1 d{k^Pn) 
d\nk 



(A34) 

for fc3 < A;*, fci » ^2 > fc* and zero otherwise. This is the 
well-known Fourier space consistency relation. 

Thus we expect the three-point correlation in real 
space to have a contribution proportional to the slope 
of the two-point correlation function. For the step po- 
tential, the slope of the two-point correlation function 
at r = 2rif diverges as Infcmax whereas S,2Tzi^^+3) di- 
verges as Infcmin- Thus, with a sufficiently small fcmin, 
we expect that the three full three-point function will be 
dominated by these squeezed contributions. Specifically, 
including the cyclic permutations which account for the 
pairings of high and low pass filtered fields we expect 



67?(r,0,0)~-21n(fc,nin??/)A^_o 



d£,2TZ 

dlnr 



(A35) 
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Bispectrum (IR) 
Consistency Rel. 
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FIG. 16. The contribution to the real space three-point 
correlation function from the infra-red terms (the second two 
terms) in Eq. ( A25 1 (dashed black curve) and the consistency 
relation of Eq. ( |A35 l (red, solid). We take Xd — 100 with an 
infrared cutoff at fcmin = O.OOI/t;/. 



near r = 2ry/ for nearly scale invariant spectra taking 
fc* comparable to but smaller than "rjj^- This explains 
the scaling of Eq. (A23). In Fig. 16 we demonstrate this 



relation numerically. Plotted is the contribution to the 
real space three-point correlation function due to the Ii 
and I2 terms of Eq. ( 48 ) and the relation in Eq. ( A35 1 . 
We take fcmin to be a factor of 10"^ smaller than the feature 
scale. The agreement is only approximate, however, this 
is possibly due to the fact that we have not omitted the 
contribution of non-squeezed triangles to ^37?. 

The In fcmin divergence is the usual IR divergence of 
a scale invariant spectrum and is cured in any physical 
observable by the finite size of Hubble radius today. The 
UV divergence from the slope of the two-point function 
at r = 2?7/ however weak, signals a breakdown of pertur- 
bation theory. In the consistency relation context, the 
assumed change in the two-point correlation in response 
to the background mode in Eq. (A27l exceeds the total 



change in the two-point correlation as a function of r. 
This is consistent with the findings of ^IVB where it 



was shown that perturbative validity places a weak con- 
straint of c/d < 10^ such that a finite step cannot have 
infinitesimal width. 



Appendix B: Slow Roll Corrections to the 
Bispectrum 

In this Appendix we show that one of the largest and 
simplest corrections to the bispectrum can be associated 
with the difference in phase between the slow roll space- 
time, and de Sitter space. This phase difference leads 
to a real component of the growing mode on super hori- 
zon scales. We show that these slow roll corrections are 
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largely responsible for the 10% error in the approxima- 
tions in 3IV| 



In deriving the expression for the bispectrum, Eq. (48 1, 



we used only the de Sitter form for the mode function 
y{x), which is related to the curvature perturbation by 



27r^ xy{x) 



(Bl) 



Here x = krj and y{x) is the solution of the Mukhanov- 
Sasaki equation, 



dx^ 



2 + g{lnx) 

1 ^ )y = 0, 



where 



/(In a;) 



/ ' 



(B2) 



(B3) 



and primes denote derivatives with respect to lux. In 
the de Sitter limit, g ~ 0, and 



y{x) = ( 1 + - I 



(B4) 



In this approximation, for modes that are outside the 
horizon, a; <C 1, the growing mode of y which yields 
the constant part of the curvature TZ is purely imagi- 
nary, while the real component decays as x^ in this limit. 
Writing the bispectrum in Eq. ( 45 1 as 



X n 



X ^ 



/ ^«^eij(ei/ -??ff)'(7^fe^7^fc^7^ 



ksl 



drj 



(B5) 



we expect that the bispectrum is dominated by the first 
term, whereas the second term should be of order ©(e^)- 
Since en ~ 0.01, one may expect that neglecting the slow 
roll corrections amounts to a small percent level correc- 
tion to the leading order bispectrum. However, as we will 
see the coefficient in front of these terms is such that the 
correction is on the order of 10% (see also [27]). 

In the limit where the inflaton is slowly rolling down 
its potential, the function ^(Ina::) is smooth, slowly vary- 
ing a nd O jus — 1). The scalar fluctuations, governed by 
Eq. (B2), begin as plane waves y oc e" in Minkowski 



space at a; 3> oo before they begin to feel the expansion 
when their wavelength becomes comparable to the hori- 
zon, a; ~ 1, where they begin to grow as y oc 1/a; for 
X <C 1. As long as g is slowly varying, the only time it is 



important is when the mode crosses the horizon. Thus 
for a given wavenumber /c, to a good approximation, we 
can Taylor expand about the time of horizon crossing, 
and to a first approximation, write 



dx"^ 



1 - 



.90 



2/ = 0, 



(B6) 



where 170 = gix = 1), is a function of k, but is time 
independent. In this limit, Eq. |B6| can be solved exactly. 
Demanding that y is asymptotically a positive frequency 
plane wave as a: — !■ cx3 implies that 



Hi'\x), 



(B7) 



where Hi'^\x) is the Hankel function of the first kind and 



(B8) 



Now, as described above in pV[ the bispectrum due to a 
step feature in the potential is dominated by modes that 
are deep inside the horizon when the inflaton crosses the 
feature. At early times, the modes are identical to the de 
Sitter modes, at leading order in x, while at late times, 
expanding the Hankel function one flnds 



xy{x) 



2 /2 

TT \ X 



. , . [2-7-ln(2)-ln(a:)] 
« + « ^ 50 

+ 0{x^gl). 



6 



(B9) 



The growing mode on superhorizon scales is no longer 
purely imaginary. There is also a correction to the am- 
plitude of the imaginary part of the growing mode, how- 
ever, note that we have already taken this correction into 



account by using the square root of Eq. ( 22 I in Eq. ( 48 ) 



which contains slow roll corrections in q 



In this ap- 
proximation, in the superhorizon limit, xy{x) is not only 
not constant but diverges as ln(x). However, it is easy to 



see that the curvature perturbation, deflned in Eq. (Bl I, 



is actually constant in this limit. The assumption that 
g = const, implies that the function / has time depen- 
dence 



f-fo 



1 



90 



In (a;) +0{gl), 



(BIO) 



to leading order in go, which is precisely the right behav- 
ior to cancel out the time dependence of Eq. ( B9 ) . 



In the limit of large momenta, corresponding to modes 
well inside the horizon, the leading order (in k and slow 
roll) correction to the bispectrum is given by 



B!^\k,,k2,k3) 



(27r)^ 

1,2 l2 1,2 

"■l "-2 ""3 



AR(fcl)AR(fc2)AR(fc3) 

4 



X T^gohiK), 



(Bll) 
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where 



(B12) 



In the same approximation as { IV we can evaluate this 
integral, to obtain 



where 



[x'^ - 9x^ + 54) sin(x) 

(-2x^ + 27x2 - 54) cos(a;) 54 



Relative to the dominant zeroth order term (the Xq term 
of Eq. ( |55[ )) the main contribution from this slow roll cor- 
rection is 7r/2 out of phase, but has the same envelope. 
Notice also that, in the limit a; — > 0, this window function 
vanishes only as X3{x) ^ 0{x), compared to 0{x'^) for 
the window functions of Eq. (55 1. This implies that, at 



sufficiently small x, this correction will eventually dom- 
inate the leading order result. This is not particularly 
surprising given that we have only considered the dom- 
inant correction in the high momentum limit. Near the 
horizon scale at feature crossing, krjf = 1, many other 
terms become equally, if not more important. 

The upper panel of Fig. 17 shows a comparison be- 



tween the analytic approximation to this first order term, 
Eqs. (Bll) and (B14| and the second term in Eq. (B5) for 
a step with height c = 10^^ and width d — 0.001. The 
prefactor for this term varies slowly over the range of 
wavenumbers considered, and we take the value Trgo/2 = 
0.1 which is near the middle of its range for this plot. 



The lower panel in Fig. 17 shows the fractional error on 
the full bispectrum when this term is added to the ana- 
lytic approximation of pV[ While the approximation is 
improved significantly at large k, the improvement gets 
steadily worse as the horizon scale is approached. This 
is to be expected, given that we have not attempted to 
calculate all of the terms that are important for modes 
that are near the horizon as the feature is crossed. 

We also note that these corrections are important for 
the power spectrum of pll| [2]. In particular, note that 
the squeezed limit consistency relation (see ^ IV C I im- 
plies that there should be a correction to the power spec- 
trum of the same order of magnitude. We leave the cal- 
culation of further slow roll corrections to future work. 

Before we end this section, we clarify one point that 
may confuse the reader. In deriving the correction 
Eq. (B9), at ffist glance it appears that the prefactor 



is simply an irrelevant global phase that will simply can- 
cel once all terms are taken into account. Indeed, the 
bispectrum is constructed out of the two-point functions 
which always involve the pair 

(7^k(ry)7^k'(^y)) =7^fc(r;)7^K77) (2^)3^3 (k-fk'). (B15) 
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FIG. 17. First order corrections to the equilateral bispectrum 
computed using the analytic approximation Eq. | |B11| | and | |B14| | 
com pare d with the numerical evaluation of the second term in 
Eq. | |B5| for a step with height c = 10~^ and width d = 0.001. 
The lower panel here shows the fractional error of the full bispec- 
trum, once this first order correction is taken into account. We 
take ■iTgo/2 = 0.1 for this plot. 



However, if one were to choose to cancel the phase in each 
propagator, one would of course obtain the same answer. 
While the growing mode would be purely imaginary out- 
side the horizon, at very early times the mode would look 
as if it initially started with a phase shift relative to the 
de Sitter space modes, leading to the same correction we 
have calculated here. 



Appendix C: Separability 

Separability of the bispectrum into products of func- 
tions that depend only on the individual fc^'s is desirable 
in that the angular bispectrum can then also be con- 
structed from separate ii calculations using the full radi- 
ation transfer function instead of the flat-sky Sachs- Wolfe 
expressions used in the main text. 

The formulation of the eigenmode decompositions of 
Fergusson, Liguori and Shellard [35] IWI which allow for 
the projection of nonseparable bispectra onto a (band- 
limited) complete basis of separable functions has largely 
ameliorated the need for the bispectra themselves to be 
separable. However, for a polynomial basis of degree n, 
one can only accurately fit a function with n zero cross- 
ings, and thus, it seems that projection of the highly 
oscillatory functions considered here might prove to be 
very inefficient with this method. Polynomials, however, 
are not the only such basis that has been proposed, the 
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oscillatory basis of Meerburg [30] would perhaps be more 
suited for this type of bispectra. Nonetheless, in this Ap- 
pendix, we point out that the analytic form of the bis- 
pectrum derived above, is to a very good approximation, 
separable. 



The bispectrum from Eqs. (48) and (54)-(55 



to be inseparable due to the fact that the damping factor 
appears to depend on the perimeter of the triangle in 
momentum space 



form. First note that at a- ^ 1 



lim V(x) = 2xe"="(l + 



(C3) 



which is separable. In the opposite limit, the damping 
factor vanishes as whereas the expansion goes as Ax. 
To fix this problem we take 



V{x) 



sinh X 



3 

i=l 



(CI) 



where Xi = kirif(nd/2^/2eQ). This cannot be exactly 
written in the form 



Vix) = Y[V,{a 



(C2) 



i=l 



We can, however, approximate the damping factor in this 



2?(x)sa2xe""^(l + e 



-4x 



(C4) 



which is of the form of a sum of separable components 
since = Y[i^~^'- The remaining inseparable 1/K" 
type factors that appear in the window functions can 
be written in separable form by introducing Schwinger 
parameters as described in |41j . However, we note that 
the terms that require this treatment are subdominant, 
and to a reasonable approximation away from K = 0, the 
bispectrum is dominated by terms that are separable. 
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